Theory of spin response in underdoped cuprates as strongly fluctuating d-wave 

superconductors 
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We study the spin dynamics in underdoped cuprates at low 
temperatures by considering them as quasi two dimensional 
d-wave superconductors with strong quantum phase fluctua- 
tions. An effective theory of spin degrees of freedom of nodal 
quasiparticles coupled to vortex defects in the phase of the 
superconducting order parameter is formulated. It represents 
the minimal extension of the QED3 theory of the pseudogap 
phase into the superconducting region. The theory predicts 
a single superconductor-spin density wave (SC-SDW) phase 
transition, without coexistence between the two phases. At 
the transition, which may be fluctuation induced first-order, 
vortices condense (and SC is lost) and the chiral symmetry 
for fermions simultaneously breaks (and SDW is formed). We 
compute the spin-spin correlation function in the fluctuating 
superconducting state, and explain the evolution of the spin 
response with energy in underdoped Y Ba2Cu30e+ x observed 
in neutron scattering experiments. In particular, we predict 
that at very low frequencies uj ~ (l/10)ov e s, with uj res be- 
ing the energy of the 'resonance' at Q = (n, 71"), (weak) spin 
response should become narrowly peaked at four diagonally 
incommensurate wave vectors that span between the nodes of 
the superconducting order parameter. These peaks represent 
the inherent collective mode of the phase fluctuating d-wave 
superconductor, the condensation of which would bring about 
the SDW order in the pseudogap phase. Our interpretation 
of the resonance suggests that it should become more elu- 
sive in the superconductors with lower T c , in accord with its' 
conspicuous absence in La2- x Sr x Cu04- 



I. INTRODUCTION 

Underdoped high temperature superconductors are 
highly anisotropic, quasi two-dimensional materials, in 
which thermal and quantum fluctuations of the phase 
of the superconducting order parameter (OP) should be 
important well below the pseudogap temperature. It has 
recently been shown that a phase fluctuating d-wave su- 
perconductor (dSC) at zero temperature (T = 0) is in- 
herently unstable towards the formation of the very weak 
spin density wave (SDW) state with loss of phase coher- 
ence [1], [2]. This follows from the observation that the 
standard Fermi liquid theory for the BCS d-wave quasi- 
particles, besides the usual spatial symmetries, also ac- 
quires an additional continuous ('chiral' [1]) symmetry at 
low energies. This symmetry becomes dynamically bro- 
ken when the superconducting phase coherence is lost 
via proliferation of vortex defects. The theory of the 
phase incoherent (pseudogap) state in this formulation 
becomes closely related to the three dimensional quan- 



tum electrodynamics (QED3) [3], in which the coupling 
constant ('charge') is proportional to the vortex conden- 
sate. In the crudest approximation, the fluctuations of 
the gauge-field may then be neglected in the supercon- 
ducting state, which leaves quasiparticles as sharp low 
energy excitations. When the vortices condense and the 
charge in the QED3 is finite, the chiral symmetry of the 
dSC becomes dynamically broken, in analogy with the 
well known phenomenon in particle physics [4]. In the 
present context this translates into weak SDW ordering, 
with confined spin- 1/2 (spinon) excitations [1]. 

Several conceptual and practical issues of direct rel- 
evance to cuprates arise naturally in such an approach 
to underdoped high temperature superconductors. Can 
the SC and the SDW long-range orders coexist? Such a 
coexistence appears rather generically in various mean- 
field treatments of microscopic models of cuprates [5], 
[6] . Several intriguing recent experiments [7] , [8] , [9] , also 
suggest the coexistence of the two inimical orderings. If 
there is a single dSC-SDW quantum phase transition, on 
the other hand, what should be its universality class? 
What are the effects of the virtual vortex fluctuations 
inside the superconducting state? These questions all re- 
quire a better understanding of the quasiparticle-vortex 
interaction, particularly in the superconducting state. 

A theory of low-energy d-wave quasiparticles coupled 
to phase fluctuations of the superconducting OP which 
puts the earlier physical ideas of [1], [3] into a consis- 
tent field-theoretical framework has recently been pro- 
posed in [10], [11]. One result of that study is that, at 
least within a certain approximation [11], in the super- 
conducting state spin of quasiparticles may be considered 
asymptotically decoupled from their charge at low ener- 
gies. Physically this is just what is expected, and may 
be understood as being a consequence of the long range 
phase coherence in a singlet superconductor [12], [13]: the 
condensate screens perfectly the charge of quasiparticles, 
effectively transforming them into neutral, spin-1/2 ex- 
citations. 

In this paper we wish to take advantage of this state 
of affairs and consider a simpler effective T = the- 
ory (Eq. (1)) for the spin sector of the fluctuating dSC 
only. We assume that in underdoped high temperature 
superconductors the amplitude of the superconducting 
OP may be considered frozen well below the pseudogap 
temperature T* , but that the quantum fluctuations of 
its phase will inevitably arise sufficiently near half filling 
when the interaction between electrons becomes strong. 
In the phase incoherent (pseudogap) state our effective 
theory reduces to the previously studied QED3 [1], [3], 
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with its concomitant SDW instability. Deep in the dSC, 
on the other hand, the theory becomes equivalent to the 
three dimensional Thirring model, with an (irrelevant) 
short range interaction between the electrically neutral 
spin- 1/2 excitations, which we will, as traditionally, call 
spinons. The effective theory in Eq. (1) connects, and 
describes the region in between, these two limits, and es- 
tablishes a physical link between the phase fluctuations 
and the magnetic excitations in undcrdoped cuprates. 

Having the low energy theory for the spin degrees of 
freedom we proceed to derive several interesting conse- 
quences from it. First, we show that the condensation 
of vortices (i. c. the loss of superconductivity) and the 
breaking of the chiral symmetry for fcrmions (the SDW 
instability) in the theory coincide. The region of coexis- 
tence between the SC and the SDW long-range orders is 
not found. Second, we find that the quantum dSC-SDW 
phase transition could be fluctuation-induced first-order, 
in agreement with our earlier study of the full theory that 
included charge [10]. Finally, and maybe most impor- 
tantly, we calculate the spin dynamics induced by vortex 
fluctuations inside the superconducting state, in the lead- 
ing approximation. The detailed evolution of the spin 
response with energy observed in the neutron scattering 
experiments on YBa 2 Cu s O &+x (YBCO) [14], [15], [16] 
is explained in terms of four diagonally incommensurate 
'mother' peaks, that should become discernible at lowest 
frequencies. These weak and narrow peaks are centered 
on the wave vectors that span between the diagonally op- 
posite nodes of the superconducting OP, and represent 
the collective particle-hole mode of the phase fluctuating 
dSC. The SDW instability obtained previously with the 
QED 3 theory can be understood qualitatively as the con- 
densation of this collective mode. The famous 'resonance' 
at the commensurate wave vector Q = (tt, tt) is here ob- 
tained as an overlap between the all four mother peaks, 
which dominates in the spin response within certain win- 
dow of energies. This energy window shrinks with the 
superconducting transition temperature T c , which may 
explain why the resonance has not been observed in 
La 2 - x Sr x Cu0 4 (LSCO) [17]. 

The paper is organized as follows. In the following sec- 
tion we introduce the effective theory of spin of quasipar- 
ticles, and discuss its physical motivation and the formal 
justification. In sec. Ill we show how our theory avoids 
the coexistence of the SC and the SDW orders. In sec. 
IV we briefly discuss the nature of the dSC-SDW quan- 
tum phase transition. Spin response is discussed in sec. 
V. The discussion of our results and the relation to other 
approaches is presented in the concluding section. 

II. THEORY OF SPINONS AND VORTICES 

Consider the quantum mechanical (T = 0) action 
for the low energy quasiparticles in the two-dimensional 
phase fluctuating dSC, S = J d 3 xC, x = (t, f), r is the 



imaginary time, and 

N 

C = Y,^il^ ll -ia ll )^i + -a-{V x A) + (1) 

|(V-zI)$| 2 + M 2 |$| 2 + ^|<I>| 4 . 

The fluctuating complex field $ describes the vortex 
loops, and ($) = implies that all such loops are of 
finite size, and thus the superconducting phase coher- 
ence. $ should be understood as being dual to the stan- 
dard superconducting OP, and ($} ^ means that vor- 
tex loops proliferated and that superconductivity is lost 
[18]. Two (N = 2) four-component Dirac fermions de- 
scribe the gapless, neutral, spin-1/2 (spinon) excitations 
near the four nodes of the superconducting order param- 
eter at ±Ki t 2, one Dirac field for each pair of diago- 
nally opposed nodes, as defined previously in [1]. a is 
the gauge-field that results from absorbing the singular 
part of the superconducting phase due to vortices into 
the spinon fields [3], [1]. 7^, jti = 0,1,2 are the Dirac 
gamma matrices [1], and {7^,7^} = 2<5 M „. Finally, A 
is an auxiliary (Chern-Simons) field which facilitates the 
statistical spinon-vortex coupling, as will be explained 
shortly. Tuning parameter fj, 2 controls the magnitude of 
vortex fluctuations, and we assume it to be related to 
doping, /i 2 oc (x — x c ), with x c being the critical doping 
in the underdoped regime. The coupling b > describes 
the short-range repulsion between the vortex loops. Wc 
work in units in which Ti = ks = 1, and we have set the 
two characteristic velocities of the nodal quasiparticles 
vp = va = 1 in (1), for simplicity. The amplitude of 
the superconducting order parameter is assumed frozen, 
and Eq. (1) is supposed to describe only the spin and 
the vortex degrees of freedom well below the pseudogap 
temperature T* , as in the previous work [1]. 

One may physically understand the form of the the- 
ory (1) as follows [19]: the exact integration over the 
gauge-field a, which appears only linearly in (1), leaves 
the spinon-vortex coupling iA ■ J$, where VXi4 = irJy, 
and J* and J$ are the spin and vortex current den- 
sities, respectively. Vortices and spinons therefore see 
each other as sources of a fictitious magnetic flux, and 
circling around a vortex with a spinon (or vice versa) 
leads to a phase change of tt. This is precisely the statis- 
tical (Aharonov-Bohm) part of the interaction between 
vortices and quasiparticles [3], [1], which couples only to 
spin. The regular part of the phase, which provides the 
standard Doppler shift of quasiparticle energies and cou- 
ples only to charge, can be 'gauged away' [13] at the price 
of introducing an additional, ultimately irrelevant inter- 
action term. Absorbing this regular part of the phase 
into quasiparticles converts them into spinons [13], rep- 
resented by ^ in Eq. (1) [1]. The final point is that at 
low energies and in the superconducting state the spin 
of quasiparticles may be considered decoupled from its 
charge [11]. Essentially, this is because one can consider 
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the superconducting condensate to carry all the charge 
in the phase coherent state. Similar conclusion may also 
be reached within the gauge theory formulation of the 
t-J model [20]. 

The separation between spin and charge allows one to 
study spin excitations separately from the charge in the 
fluctuating dSC at low energies. Eq. (1) may there- 
fore be understood as an effective theory for the spin of 
quasiparticles. One can also derive it from the full the- 
ory which includes charge [10], [11] by simply omitting 
the additional Chern-Simons field for the charge degrees 
of freedom. 

Finally, in writing (1) the anisotropy between the 
Fermi and the OP related velocities, vp 3> va, and 
all the local interactions between quasiparticles [2] have 
been neglected, as irrelevant at low energies [21], [22]. 
We will, however, find it necesary to restore the velocity 
anisotropy later when we discuss the spin response, where 
this feature will turn out to play an important role. 

Although in this work we are concerned with the vor- 
tex fluctuations in the superconducting state, which is a 
complementary problem to the one studied in the previ- 
ous papers on the QED3 theory of the pseudogap phase 
[1], [3], let us nevertheless mention how the theory (1) 
is related to the QED3 in the mean-field approximation 
for the vortex field. Entirely neglecting the fluctuations 
in $, in the superconducting phase (($) = 0), the inte- 
gration over A would constrain V x a = 0, and spinons 
become free. When vortices condense (($) ^ 0), on the 
other hand, the gauge-field A acquires a mass via Higgs 
mechanism, |($)| 2 A 2 . The Gaussian integration over A 
produces then the Maxwell term for a, ~(Vx a) 2 /|($}| 2 . 
Together with the Dirac Lagrangian this constitutes the 
QED3 for spinons, in which the fermions' chiral sym- 
metry generated by 73 and 75 is dynamically broken by 
the generation of the small mass term ~ M^^i [4]. In 
the present context such a mass M ~ (^i^i) ~ K^)! 2 is 
proportional to the weak SDW OP, with the (incommen- 
surate) ordering wave vectors ±2Ki, 2 that connect the 
diagonally opposed nodes of the superconducting OP [1] . 

Our aim will therefore be to understand how the the- 
ory (1) interpolates between the limits of QED3 in the 
pseudogap, and the non-interacting theory of spinons in 
the superconducting phases. In particular, we want to 
understand the evolution of spin dynamics with progres- 
sive quantum disordering of the dSC, which is assumed 
to correspond to undcrdoping. We begin with the issue of 
possible coexistence of the SDW and the dSC orderings. 



III. COINCIDENCE OF THE SC AND THE SDW 
TRANSITIONS 

In the mean-field approximation, the vortex condensa- 
tion and the SDW transition in the theory (1) coincide. 
It is not obvious that this feature survives the inclusion 
of the fluctuations, and it is conceivable that the SDW 
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FIG. 1. Schematic behavior of the inverse gauge-field prop- 
agator 1/Gaa(p, m) as a function of the momentum in the SC 
(m 2 > 0), critical (m = 0), and the pseudogap (m 2 < 0) 
phases of the theory (1), from top to bottom. At large mo- 
menta the behavior is always linear. At low momentum, how- 
ever, the statistical gauge field is massive (~ m) in the SC 
phase, and has a Maxwell term (~ p 2 ) in the pseudogap 
(incoherent) phase. 



transition may occur within the superconducting phase 
[1] , [5] , [6] , [23] . To examine this issue we first ask what 
would be the exact propagator for the gauge field a if 
there were no fermions, i. e. when N = in (1). We then 
approximate the interacting action for a that would re- 
sult from the integrations over $ and A with the effective 
Gaussian term that reproduces that exact propagator. In 
the Landau gauge such a propagator is 



n a , \ 7r 2 n AA (p,m) 



(2) 



where Haa(p, tn)(&nv — p^Pv) is the transverse current- 
current correlation function in the |$| 4 theory. Here 
m 2 = /1 2 + 0(b) is the fully renormalized "mass" of the 
vortex field. In general, Haa(p, m) = npF + (m/p), for 
m 2 > (superconductor). To the lowest order in b 



F+(z) = ---[(4z 2 + l)arctan(^) - 2z] + 0(b). 



(3) 



Here we generalized our model to the one with n complex 
vortex fields, n = 1 being the case of physical interest. 

We are now in position to study the chiral symmetry 
breaking for fermions with (2) serving as the bare (with- 
out fermion polarization) gauge-field propagator in the 
QED 3 . Consider the standard large-N Dyson equation 
[4] for the fermion self-energy 



4 7 (27T 



d 3 p_ Gga^ v (p - q,m)Yl(p) 



P 2 + S(p) 2 



7, (4) 



Here [G aa ^ v (p, m)} 1 = U^ a ^(p) + [G° aa ^(p,m)} \ 
with the one-loop fermion polarization H^ a ™(p) = 
(Np/8)(S^ — Pfj.Pi,)- It is useful to consider first the 
point of the superconducting phase transition m = 0, 
where Eq. (3) implies U AA (p,0) = np(l/16 + 0(b c )), 
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where b c is the fixed point value of b in d — 3 [24] . Linear 
dependence of Haa{p, 0) on momentum is an exact result 
[25] . Inserting this into the Dyson equation we find that 
at m = the effect of vortices is only to increase the co- 
efficient N — ► N e ff = N + 128/(7r 2 n) + 0(b c ). Recalling 
that the non-trivial solution of the Dyson equation (4) 
exists only for N eff < N c = 32/ tt 2 [4], for N = 2 we find 
that chiral symmetry at m = would be already broken 
only for n > n c = 10.44(1 + 32/(97r 2 n)), where we have 
also included the known 0(b c ) correction to Haa(p, 0) in 
the large-n approximation [24]. Since n = 1 in the phys- 
ical case we conclude that right at the superconducting 
critical point SDW order is most likely absent. 

Little further thought shows that, if correct, the above 
result implies that the SDW order is absent for all 
to 2 > 0. Indeed, this is to be expected, since G° a (p, to) < 
G° aa {p, 0), and a is only stiff er in the superconducting 
phase. To prove this, assume that for m 2 > a non- 
trivial solution of the Eq. (4), £(<?), does exist. Such a 
solution would then satisfy 



t(q) < \tt 



(p-g,0)£(p) _ 



(27T) 



P 2 + S(p) 2 



On the other hand, we already established that there is 
only a trivial solution of the Dyson equation at to = 0. 
This means that assuming a candidate function E(p) and 
inserting it under the integral in the Dyson equation (4) 
for to = will produce only a smaller new £(</), since un- 
der iterations the physically acceptable self-energy must 
approach the trivial solution. This further implies that 
any £(<?) from the domain of attraction of the trivial so- 
lution at to = has to satisfy the inequality opposite to 
(5). A non-trivial £(g) at to 2 > therefore can not exist 
if it did not exist at to 2 = 0. 

For to 2 < superconductivity is lost, and T1aa(p, to) = 
npF_(\{&)\ 2 /p), with F_(z) -> 1/16 for z < 1, and 
F-(z) = 2z, for z 3> 1 (Higgs mechanism), to the lead- 
ing order. The behavior of the gauge field propagator 
is depicted in Fig. 1 . Since G° aa (p, to) > G° aa (p, 0) for 
to 2 < 0, the above argument no longer applies. In fact, 
it is easy to see that there is immediately a non-trivial so- 
lution of the Eq. (4) when to 2 < 0. Since £(<?) > only 
for q <~ |(3>)| 2 and smaller, |($)| 2 serves as the effective 
ultraviolet cutoff in the Eq. (4). The Dyson equation 
then reduces to the standard one in the QED 3 [4]. For 
N = 2 < N c , T,(q = 0) ~ |($)| 2 , and the chiral sym- 
metry is dynamically broken. There is no intermediate 
(quantum disordered) phase in between the SDW and 
the dSC in the theory (1), unless the exact value of N c 
in the QED3 is actually less than two [26], [27]. 



IV. NATURE OF THE TRANSITION 

The above argument strongly suggests that there is a 
single dSC-SDW transition in the theory (1), but does 
not say anything about its nature. This issue in the 



full theory with both spin and charge was addressed in 
[10], and here we will provide a brief alternative cal- 
culation in our effective theory that supports the ear- 
lier results. To study the phase transition it is better 
to proceed in the opposite direction, and integrate over 
the fermions first. This gives dynamics to a via the 
fermion polarization bubble. If we neglect the quartic 
and the higher order terms and perform the Gaussian 
integral over a, the result is the Maxwell-like term for 
A: (A^p)A u {-p)) = {NTr 2 /8\p\){5^- Pl , Pv ). Note that 
since the inverse of the above average is non-analytic at 
p = 0, it can not renormalize [28], and therefore the 
number of fermions TV represents an exactly marginal 
coupling. Assuming a constant $ in (1) we may further 
integrate A to find the energy per unit volume to be 



(6) 



where we have rescalcd $ 2 , with the ultraviolet cuttof 
A implicit in (1). Standard analysis of the Eq. (6) 
shows t hat t here is a discontinuous transition for N > 
(4/7r)^/26/A). Using the lowest order fixed point value 
b/A = (2tt 2 /5)(4 -d) + 0((4 - df ) in d = 3 as a crude 
estimate of this bound, we find the first-order transition 
for N > 3.58, in rough agreement with [10]. Of course, 
this conclusion is to be trusted only for N 3> 1, when 
the first-order transition occurs at a large /x 2 , at which 
our neglect of fluctuations in $ in arriving at the Eq. (6) 
becomes justified. For JV < 1, on the other hand, the 
phase transition in the theory (1) should be continuous, 
with weakly modified XY exponents: v = v xy + O(N), 
V = Vxy + 0(N) [10], [19]. The situation is reminiscent of 
the Ginzburg-Landau superconductor [25], with N ^> 1 
analogous to the strongly type-I, and N <C 1 to the ex- 
treme type-II case. Of course, in dSC N — 2 and fixed, 
so the issue is whether this falls into the type-I or into the 
type-II regime. In this context it may bc interesting to 
note that T c does appear to drop discontinuously to zero 
at the critical doping in LSCO [17], and very steeply, if 
not even discontinuously, in YBCO [29]. 



V. SPIN RESPONSE IN THE 
SUPERCONDUCTOR 

A. Derivation of the response function 

Finally, we turn to spin dynamics deep inside the dSC, 
when m/b 3> 1. For small momenta, p <C Aj-f, w 5m, we 
find that G° aatflv (p,m) = ((n/24m) + 0(p 2 ))(6^-p^p v ). 
To keep the algebra simple, in the following wc will ap- 
proximate the gauge-field propagator with the constant 
mass term, corresponding to the leading contribution in 
the above. The retention of the full form of G° aa (p, to) 
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should not qualitatively change the results, although it 
is more than likely that there will be some quantitative 
differences. 

Integrating out such a massive a in the superconduct- 
ing phase leads to the effective 2 + 1 dimensional Thirring 
model [30] for spinons 



72m 



(7) 



with the summation over repeated indices and the ul- 
traviolet cutoff A.Th assumed. Using the Hubbard- 
Stratonovich transformation [31] we can rewrite this ex- 
actly as 

C = Va^i + —Tr[Mj>M»] + TrlM^^^i] (8) 



0.06 



0.02 



-0.02 



-0.06 




-0.06 -0.04 



. 04 0.06 



Making an ansatz M l J (x) — (M(x)/3)Sijj l _ l at the saddle 
point, and then integrating out the fermions yields 

/ 94777 
d 3 x{— M 2 {x) - Tr[ln( 7 ^ - M{x)))]. (9) 

Expanding further in powers of M(x) we finally write 



(10) 



where M 2 (q) = M(q)M(-q). 

Using the definition of the Dirac fields in terms of the 
electron creation and annihilation operators [1] it readily 
follows that 

2 



i=l 



where the vectors ±K i7 i = 1,2 denote the positions of 
the four nodes of the d-wave order parameter, as before. 
(M(q = 0)} is therefore the static SDW OP [1], which, 
of course, was the motivation behind our saddle-point 
ansatz. Similarly, 



(M 2 (q,co))- 



48iVm 



= (12) 



i=l 



(S z (q- 2Ki,tu)S z (-q + 2Kt, -co))). 



Both Eqs. (11) and (12) represent exact equalities. 

Since the dSC is rotationally invariant, the spin- 
spin correlation function is diagonal in indices a, /3, 
(S a (k,u))S0(—k,—u))) = x(k,co)6 a p. We may therefore 
finally deduce that the imaginary part of the spin re- 
sponse function = x'{k,uu) + ix"(k,uj) in the 
phase fluctuating dSC is 



FIG. 2. Contour plot of the x" a t a l° w frequency 
uj — 6meV, where the spin response peaks at four diagonally 
incommensurate 'mother' peaks. The origin corresponds to 
(it, it), and numbers are in units of 2tt. For other parameters 
see the text. Brighter shade denotes a higher intensity, y" 
vanishes in the black regions. 

_ 1 O AT 777 

X "(2K i ± q,u) = ( flm(M 2 (q,w)). (13) 

Neglecting the quartic and higher order terms in Eq. 
(10) (which may be formally justified by assuming N 3> 
1), and analytically continuing to real frequencies iuo — > 
u, in the Gaussian approximation we finally obtain 



X "(2K i ± q, uj) = — 2 \ q ' V 2 14 

7T (64to/vt)^ + or — q z 

where we have written the cutoff At/i ~ 5m as Ath — 
(16/7r)m, for convenience. Eq. (14) is our central result. 
Its' most important feature is that the relativistic form of 
the quasiparticle dispersion near the nodes forces x" to 
vanish below a certain momentum-dependent frequency. 

First, let us check that deep in the superconducting 
phase we recover the result for the free quasiparticles. 
Taking m — > oo in the Eq. (14) 

N 



lim X "(2X, ± q,u) = -Q(co 2 - q 2 )y/oj^f, (15) 

m— >oo 54 

in agreement with [32] and [33] for N = 2. For q = 0, 
for example, x" ~ M> an d the spin response at ±2Ki is 
completely suppressed, as expected. 

Next, we connect the uniform vortex susceptibility m 
to the measurable superconducting T c . Assuming a con- 
tinuous dSC-SDW transition, m 2 oc (^i 2 — y(i 2 ) 7 , with 
7 = i/(2 — 77) characterizing the quantum critical point; 
also, T c cx ([i 2 — ii 2 ) zp '. Since the dynamical critical ex- 
ponent in the theory (1) is z = 1, m oc T c 1_ ^^ 2 \ Suffi- 
ciently away from the critical point we may neglect the 
anomalous dimension 77 and find the result dictated by 
the (engineering) dimensional analysis, m = c(-kT c / d>4) , 
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where c is a number. This would remain true even if the 
transition is weakly first order. 

On the other hand, retaining the charge degrees of 
freedom one finds that the T = supcrfluid density is 
p s f w 24to/-7t, in the underdoped regime [11]. This al- 
lows us to estimate the number c by utilizing the Ucmura 
scaling [34], by which in our units p s f rts 2.1T C in under- 
doped cuprates [35]. This number is known not to be 
truly universal, but due to structural similarities between 
the cuprates the variations from one material to another 
are reasonably small [35]. The results that follow will 
at any rate turn out to be quite insensitive to its exact 
value. This finally yields c w 5.6, i. e. 



This should be understood as only a crude estimate; 
changing the cutoff Ax/j, for example, will alter the nu- 
merical value of the constant c. This inherent limitation 
of the present calculation underlines the need for a better 
treatment of the gauge field term in the superconduct- 
ing phase, that would improve upon our constant mass 
approximation. This would, however, make our calcula- 
tion significantly more technical, and obscure somewhat 
the simple physics we wish to address. Keeping these 
cautionary remarks in mind, we proceed to observe the 
consequences of Eq. (14) in the following section. 

B. Evolution with frequency 

We look for the maxima of the x"(<f, w) in the mo- 
mentum space at a fixed u>. There are four regimes to 
discern: 

1) For low frequencies u>, the maximum values of \" 
are located at four "diagonally" incommensurate wave 
vectors ±.2K\^ (i. e. at q = 0). In Fig. 2 we plot 
such four well separated peaks at 2K\, —2K\ + (2tt, 2tt), 
2K 2 + (2n, 0), and — 2K 2 + (0, 2ir), where we have used the 
periodicity provided by the underlying lattice to shift all 
the peaks into the vicinity of the commensurate wave vec- 
tor Q = (7r,7r). For illustration, we assumed T c = 67 K 
[16], vf/va — 7 [36], vf = 1.2eVA [37], and estimated 
2Ki = 0.465(27r, 2ir) from [14]. These 'mother' peaks are 
narrow and of low intensity, since both their width and 
the magnitude are controlled by the frequency. 

2) As the frequency increases both the the intensity 
and the width of the mother peaks in the momentum 
space grow. Since vp 3> v\ they first start to overlap 
at four 'parallel' positions (Fig. 3). This way the su- 
perposition of the mother peaks leads to appearance of 
the maximum response at four incommensurate 'paral- 
lel' wave vectors. We emphasize that this type of in- 
commensuration has nothing to do here with any sort of 
one dimensional ('stripe') ordering, as often assumed in 
literature. 
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FIG. 3. The same as Fig. 2, but at a larger frequency, 
u) — 15meV . Note that the maximum of the spin response 
is now found at the 'parallel' incommensurate wave vectors, 
because of the overlap of the mother peaks. 

3) With a further increase of frequency the four initial 
peaks start overlapping at Q — (tt, tt), and for a while the 
'commensurate' response dominates (Fig. 4). With the 
above parameters we find the response at (tt, tt) starting 
developing at ~ GOmeV (see Fig. 6). The energy of 
the 'resonance' at which the response at Q = (tt, tt) is 
maximal is oj res « 68meV. The overlap of all four mother 
peaks makes the response large in the whole interior of 
the square in Fig. 2. 

4) Finally, at the largest frequencies the maximum in 
Eq. (14) shifts to a \q\ ^ 0, which implies a weak redistri- 
bution of the commensurate peak to four 'parallel' incom- 
mensurate positions again (Fig. 5). However, whereas 
the position of the 'parallel' incommensurate peaks at 
Fig. 3 was independent of frequency at low frequencies, 
now it increases with frequency. 

C. Comparison with experiment 

Let us review some of the features of neutron scat- 
tering experiments on underdoped cuprates pertaining 
to our discussion. We will concentrate on the under- 
doped YBCO at low temperatures, where the most de- 
tailed study of spin response as a function of frequency 
and wave vector exists [16]. 

1) At low frequencies (w < 40meV) the maximum re- 
sponse is at four 'parallel' incommensurate wave vectors 
[8], [16]. The location of these peaks is very weakly de- 
pendent, or almost independent, on frequency [16]. 

2) At a frequency uj <~ AOmeV the maximum shifts to 
the commensurate position at (tt, tt) , and peaks there at 
a slightly larger energy. This is the so-called 'resonance'. 

3) At a higher frequency ui ~ 50meV the resonance 
becomes flat, and two weak 'parallel' incommensurate 
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FIG. 4. Spin response at w = 72meV, near the resonance 
energy of u) res — 68meV. The overlap of all four mother 
peaks from Fig. 2 has now produced the largest response at 
the origin, i. e. at the commensurate wave vector Q = (tt, tt). 

peaks can be discerned again. This time, however, their 
incommensuration increases with frequency. 

4) The energy of the resonance decreases with the su- 
perconducting T c [14], [15]. 

5) At Q = (tt,tt) there is a 'spin gap': the commen- 
surate response is strongly suppressed below a certain 
energy oj sg . At high frequencies the response also contin- 
uously decreases, leaving a maximum at the 'energy of 
the resonance'. 

6) Overall spin response increases with frequency at 
low frequencies, peaks around the energy of the reso- 
nance, and then decreases with further increase in fre- 
quency 

We see that the features l)-3) are at least qualitatively 
reproduced by our Eq. (14), once the lattice periodicity 
and the velocity anisotropy are taken into account. Also, 
if we assume that Vf and va are only weakly doping de- 
pendent, the only energy scale left is the gauge field mass 
to, which is connected to the doping dependent super- 
conducting T c . The observed scaling of resonance energy 
with T c in underdoped samples therefore automatically 
follows from our theory. Note that this feature otherwise 
may appear quite counterintuitive: if one believes that 
in underdoped cuprates T c is the temperature where the 
phase coherence sets in, while the higher and increasing 
pseudogap temperature T* is related to the amplitude 
of the OP, the observed connection between the spin re- 
sponse and the superconducting phase coherence would 
seem rather mysterious. 

The resonance at Q — (tt,tt), in our notation corre- 
sponds to some finite qo = Q — 2K\. This implies that 
as a function of frequency, at T = \" vanishes be- 
low the spin gap u sg w vpqo ~ 60meV, which is in 
our theory essentially determined by the incommensura- 
tion of the mother peaks and the magnitude of vf (Fig. 
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FIG. 5. Spin response at a high frequency ui = HOraeV. 
Note the weak redistribution of the intensity from the 'reso- 
nance' in the center to the 'parallel' incommensurate peaks, 
which now will also acquire an upward dispersion. 

6). While our estimate appears to be too large, our re- 
sult naturally accounts for the existence of the spin gap. 
Also, the maximum response in Eq. (14) is at the energy 

uj res = ^JiUgg + (cT c ) 2 , which decreases with decreasing 

T c , as observed [14], [15]. The experimental data of Dai 
et al. [14] are well reproduced by this simple formula with 
a just slightly smaller constant c than in Eq. (16), and 
lu scj = 22meV (Fig. 7). As T c — > 0, therefore, the com- 
mensurate peak energy u; res — > uo sg ^ 0. We thus predict 
that the spin gap, and the T c — ► extrapolated energy 
of the resonance are essentially the same. Remarkably, 
the spin gap that can be estimated from a different set 
of data on YBCO (Fig. 17 in [15] ) is quite similar, 
uj sg w 20meV, consistent with our interpretation. At 
large energies, \" should behave as ~ 1/w (see Fig. 6), 
which also appears to be in a qualitative, but not quite 
a quantitative, accord with the data. 

Finally, we predict a new feature of the spin response at 
the lowest frequencies: four weak and narrow 'diagonally' 
incommensurate mother peaks, the energy of which van- 
ishes with T c . Their low intensity and sharpness make 
them rather elusive, and indeed they have yet to be ob- 
served. Nevertheless, the fact that the evolution of the 
spin response with frequency as observed in underdoped 
YBCO can be straightforwardly understood in terms of 
these makes one optimistic about their existence. Wc 
find that the mother peaks should resolve as in Fig. 2 
at a frequency ~ u) res /10. Our optimism is further sup- 
ported by the fine details of the measured spin response, 
which appear so far to be in accord with our picture. For 
example, Eq. (14) predicts that the response function at 
the 'parallel' peaks must always be less than twice larger 
than that at ±21?!^ at the same frequency. Closer in- 
spection of the data in [14] or [16], for example, conforms 
to this expectation. 
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FIG. 6. Spin response function \" ( m arbitrary units) at 
the commensurate wave vector Q = (ir, n) as a function of 
energy (in units of vf = 1.2eVA). Note that the spin response 
vanishes below certain frequency ('spin gap'), sharply peaks 
just above it, and behaves like ~ 1/lu at large frequencies. 

VI. CONCLUSION AND DISCUSSION 

We have proposed an effective theory for the low- 
energy spin excitations in the strongly fluctuating d-wave 
superconductor. The theory extends the QED 3 of the 
pseudogap phase into the superconducting phase. It sug- 
gests that there is a single dSC-SDW transition, which 
may be fluctuation induced first order. We derived ap- 
proximately the spin response function of the fluctuating 
d-wave superconductor, and used it to qualitatively ex- 
plain the neutron scattering experiments in underdopcd 
YBCO. 

Besides standing in clear opposition to most of the 
mean-field theories of the high temperature superconduc- 
tivity [5] , [6] , our conclusion that there is no coexistence 
region between the dSC and the SDW is also in contrast 
to the result of Pereg-Barnea and Franz [23] , who studied 
a closely related theory to (1). In their model spinons are 
coupled to the massive gauge field in the superconducting 
phase, so that after integration over the gauge field the 
model would become similar to the Thirring model in our 
Eq. (7) . This leads to the conclusion that the chiral sym- 
metry for fermions becomes broken at a finite mass of the 
gauge field m, which implies coexistence of the SDW or- 
der and the phase coherence. One can readily understand 
this result by considering the dimensionless coupling con- 
stant in the Thirring model g — Axh/m. Indeed, if the 
cutoff Ath is kept constant, g — * g c ~ 1 at a finite to, 
and the model exhibits its well-known chiral instability. 
This result, however, in the more complete theory (1) 
is avoided by remembering that in the Thirring model, 
or in an equivalent description of [23], the cuttof needs 
to be Axh ~ to as well, because for momenta p 3> m 
the inverse gauge field propagator becomes a linear func- 
tion of p (Fig. 1). With this modification, as to — > 0, 
g — ► constant, with the constant apparently too small 
for the chiral instability to occur before the critical point 
is reached. 

The result that superconducting and antiferromag- 
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FIG. 7. The energy of the resonance (in meV) vs. the 
superconducting T c (in K). Data are from the Fig. 29 in 
second ref. 14, and the fit is to cj res = */ (cT c ) 2 + ui'ig (see 
the text), with c = 4.75, and u> 3g — 22meV. 

netic orders are not comfortably coexisting in under- 
doped cuprates is also in agreement with recent exper- 
iments on heat transport and c-axis penetration depth 
in clean underdoped YBCO [38], [39]. Both measure- 
ments find that d-wave nodes survive in the underdopcd 
regime, without the quasiparticles becoming gaped. If 
there would be a sufficiently strong coexisting antiferro- 
magnetic (Q = (tt, tt)) or SDW (Q ^ (it, it)) order, on the 
other hand, quasiparticles would acquire precisely such a 
gap at the nodes. It should be said, however, that the 
coexistence of SDW and SC orders is not in principle 
impossible, and can be induced by a strong favourable 
interaction term, for example, omitted in (1). Our point 
is only that vortex condensation as the dynamical mech- 
anism behind the appearance of the SDW order in the 
theory (1) does not lead to such coexistence. 

Our picture of spin response in underdoped cuprates is 
similar in spirit to that of Rantner and Wen [32] . There 
are important conceptual and formal differences, how- 
ever. Rantner and Wen considered the pseudogap regime, 
above the superconducting T c , and proposed that spinons 
may there effectively be described as coupled to a mass- 
less gauge field. Neglecting the chiral instability such a 
theory is almost universally believed to have, they pro- 
posed the algebraic spin liquid as the pseudogap phase, 
in which in the large-N approximation the spin response 
is in fact suppressed at low frequencies. Only if the num- 
ber of fermion components N is taken sufficiently small 
(smaller in fact than the N c below which the chiral in- 
stability occurs in the same calculation) does one recover 
the enhancement of the spin response. To account phe- 
nomenologically for the phase coherence in the supercon- 
ducting state Rantner and Wen added a mass to their 
gauge field, which then further suppresses the spin re- 
sponse at energies below the mass. We see that a qual- 
itatively similar mechanism operates in our theory, but 
at T = 0: the gauge field coupled to spinons is mass- 
less outside the SC, and massive inside. The form of our 
spin response function x", however, is quite different, 
which, for example, leads to reappearance of the incom- 
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mensurate peaks at high energies (Fig. 5), among other 
things. Most importantly, the mass of the gauge field 
m in the superconducting phase that was postulated in 
[32] on phenomenological grounds here becomes a result 
within the effective theory (f). 

With the parameters we used the energy of the reso- 
nance LOres ~ 70meV , not far from the generic ~ 40meV 
see in YBCO, and more recently in Bi 2 Sr2CaCu20s+s 
(BSCCO) [40], [41]. Given the crudeness of our approxi- 
mation this is an encouraging agreement. However, there 
is a considerable uncertainty in the parameters we used 
which may easily make this agreement fortuitous. For 
example, assuming vp — l.SeVA [42] would push the 
resonance energy in our calculation to ~ lOOmeV. Also, 
we do not find the (5-function response at the commensu- 
rate wave vector, but only al/w decay on the high energy 
side (Fig. 6). This is contrast to the mean-field calcula- 
tion within the t — J model of Brinckmann and Lee [6] , 
who found that the next-nearest-neighbor hopping may 
lead to a ^-function term in x" ■ Our result, however, 
seems to agree with the measurements of Fong et al. [15] 
on underdoped YBCO, where the width of the resonance 
of ~ lOmeV was observed (see Fig. 7 in [15]). It is inter- 
esting that Brinckmann and Lee, just like us, also tend 
to overestimate the energy of the resonance, and found it 
necessary to reduce the coupling J by hand by a factor 
of three in their calculation to match the experiment. 

We find the commensurate response to dominate only 
within a certain window of energies around the reso- 
nant energy, while both at lower (u < ui sg ) and higher 
(lu > l.5u> sg ) energies the peaks are located at the par- 
allel incommensurate wave vectors. In our calculation 
the resonance exists for 60meV < u < 90meV, with the 
parameters we used. For a lower T c the window of the 
dominant response at Q = (ir, ir) shrinks. This may pro- 
vide an explanation why the resonance was never found 
in LSCO, but has been seen in both YBCO and BSCCO 
which both have a higher T c , without other parameters 
being too different. Alternatively, the physics of under- 
doped LSCO may be qualitatively different, as some mea- 
surements would suggest [17]. 

We end by repeating the central prediction of this 
work: the appearance of four weak and narrow 'diago- 
nal' incommensurate peaks, the energy of which should 
extrapolate to zero with vanishing superconducting T c . 
It is the interplay of these fundamental peaks that leads 
to evolution of spin response with frequency in our calcu- 
lation in qualitative agreement with experiment. While 
some other theories [6], [32] would also lead to such 'diag- 
onal' incommensuration at low frequencies, we also pre- 
dict that the energy of these peaks should vanish with T Cl 
since they represent the intrinsic soft mode of the phase 
fluctuating d-wave superconductor. Detection of this fea- 
ture in future neutron scattering experiments would pro- 
vide strong support in favor of the theory of the under- 
doped cuprates advocated here. 
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